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SUMMARY 

We formulate the problem of finding the probability that the determinant of a matrix undergoes the 
least change upon perturbation of one of its elements, provided that most or all of the elements of the 
matrix are chosen at random and that the randomly chosen elements have a fixed probability of being 
non-zero. Also, we show that the procedure for finding the probability that the determinant undergoes 
the least change depends on whether the random variables for the matrix elements are continuous or 
discrete. Preprint submitted to Numer. Linear Algebra Appl. on June 28, 2007 
Prepared using niaauth.cis 
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1. Introduction 

Consider the problem of finding the probability that the determinant of a matrix undergoes 
the least possible change upon perturbation of one of its elements. In this paper, we present 
non-asymptotic solutions for that class of matrices and others, and we discuss the connection 
between the case where the (real) random variables for the matrix elements are continuous 
and the case where they are discrete. Denote the determinant of any matrix A by det A. For 
an (rt + 1) X (n + 1) matrix Mn+i, the expansion of det Mn+i via row i is 

det Mn+i = muMii + h rriijMij H h mm+iMm+i, (1) 

where rriij is the element of Mn+i at the intersection of row i and column j, and My is 
the cofactor of rriij. Since rriijMij is the only term in that expansion of detM„+i in which 
element my occurs, the probability that perturbation of rriij has the least effect on the 
value of detM„+i is equal to the probability that |My | is as small as possible. The cofactor 
Mij of element rriij in Mn+i is (— 1)*+^ det 5'„, where Sn is the n x n submatrix of M„_|_i 
which is obtained by deleting row i and column j. Thus |My | = | det Sn\, so the probabihty 
that detM„^i undergoes the least change upon perturbation of element rriij is equal to the 
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probability that | det 5'„| is as small as possible. In particular, perturbation of rriij has no effect 



at all on det M„+i if det Sn 
matrices Sn- 



0. In this paper, we treat the following three classes of n x n 



(i) 



(ii) 



Matrices An in which all the elements are values of mutually independent random 
variables each of which has probability r of being non-zero and probability 1 — r of 
being 0, where < r < 1. 

Matrices B„ in which all but one of the diagonal elements are set to 1 (i.e., bu = 1 
for i ^ 1, where bu is the special diagonal element), and bu and all the off-diagonal 
elements are as in (i). Thus 





1 bu bi2 bi3 


••• bin \ 




&21 1 &23 


■■■ b2n 


Bn 


bsi &32 1 


■■■ ban 




\ bnl bn2 bn3 


... 1 ) 


r which all the diagonal elements 


are set to 1 


as in (i). Thus 








/ 1 Ci2 Ci3 


• • • Cin \ 
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■■■ C2n 




C3I C32 1 


■•■ C3„ 




V Cnl Cn2 C„3 


... 1 ) 



The matrices 5„ described in (i), (ii), and (iii) above will be said to be matrices of type A„, B„, 
and C„, respectively. Every randomly chosen element of a matrix of any of these typos (i.e., 
every element of a matrix of type A„, every off-diagonal element of a matrix of type B„ or C„, 
and the special diagonal element 611 of a matrix of type B„) will be called a variable element; 
the remaining elements (those which are not chosen at random) will be called fixed elements. 
Let jC be a set such that € X, let X be a random variable whose values are in X, and let P 
be the probability function for X defined by 



P{X = 0) = 1 - r 
P{X^O) = r 



(2) 



If Z? C X and X is a discrete set of numbers, then P{X G D) can be expressed in terms of a 
probability function p as 

P{x eD) = J2 P(^). (3) 

where 

0^p{x)^l, p(0) = l-r, ^ p{x)=r 

xex-{0} 

If X is a continuous set of real numbers, then P{X G D) can be expressed in terms of a density 
function q as 

P{X €£>)=/ q{x)dx, (4) 
Jd 
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where 

^ g(a;) ^ 1, / q{x)dx = l — r, / q{x)dx = r 
J{0} Jx-{0} 

In what follows, q is assumed to be continuous on X — {0}. Then / q(x)dx — D — {x} 

Jd 

for some x G X — {0}, which implies that 

Proposition 1.1. Let X be a continuous set of real numbers such that G X, and let M„_|_i 

be an (n + 1) x (n + 1) matrix all of whose variable elements are in X. Furthermore, let Sn 
be the n x n submatrix of Mn+i which is obtained by deleting row i and column j. Then the 
problem of finding the probability that detM„+i undergoes the least change upon perturbation 
of element nriij ( of Mn+i) is equivalent to the problem of finding the probability that det Sn = 
if S„ is of type or B|^j and to the problem of finding the probability that d^itS^i — 1 if Sn 
is of type C„. 

Proof Expand the determinant of Sn via the permutation group S„ on the set {1, 2, 3, ... , n}: 

det5„= sgn{a)s„(^i)iSa(^2)2- ■ ■ Sa{n)n, 

where sgn^a) is the sign of permutation a (sgn{(T) = 1 if cr is an even permutation, and 
sgn{a) = — 1 if (J is odd). First, consider a matrix Sfi of type A:^ or B„. Claim 1 Let xq be a 
non-zero number. If Sn is of type A„ or B„, then the probability that detSn = xq is 0. Proof 
of Claim 1 If detS'„ = xq, then at least one term in the expansion of dctS'„ is non-zero, and 
every element of Sn that occurs in at least one non-zero term in that expansion is non-zero. 
Since Sn is of type A„ or B„, every non-zero term in the expansion of det 5„ contains at least 
one variable element of Sn- Because dctS'„ = .tq, the sum of all the non-zero terms must be 
xo, so the value yo of the "last" non-zero variable element is uniquely determined from all the 
others. For example, if n = 2 and 52 is of type A2, then all the elements of 52 are variable. 
Because 52 is a 2 x 2 matrix with a non-zero determinant and has at least one non-zero term 
in the expansion of its determinant, at least one of the products S11S22, S21S12 is non-zero. 
Suppose that S11S22 is non-zero. If we consider S22 to be the last element, then 

, , c- . f \f \ f \f \ . a;o + (S2i)(si2) 

detA2=xo (sii)(s22) - (S2i)(si2) = a;o =^ S22 = Vo = 

If yo ^ X, it is obvious that P{X ~ y^) ~ 0, so assume that yo £ X. Since X is a continuous 
set, it follows from (5) that P{X = yo)=0. From this result, together with the usual rules for 
computing the probability that a finite sum of finite products of numbers (such as the expansion 
of det Sn via 6„) has a certain value, it follows that the probability that dot 5„ = xq is 0. End 
of Proof of Claim 1 By Claim 1, the probability that detM„_|_i will undergo the least change 
upon perturbation of element ruij is equal to the probability that det 5„ = 0. Claim 2 If 5„ is 
of type A„ or B„ , then the probability that det Sn = and at least one term in the expansion 
of det Sn via Sn is non-zero is 0. Proof of Claim 2 If at least one term in the expansion of 
det Sn is non-zero, then every element of 5„ that occurs in at least one non-zero term in that 
expansion is non-zero. Since 5„ is of type A„ or B„, every non-zero term in the expansion 
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of dot Sn contains at least one variable element of Sn- Because det S'„ = 0, the sum of all the 
non-zero terms must be 0, so the value yo of the "last" non-zero variable element is uniquely 
determined from all the others. Just as in the proof of Claim 1, the probability that that last 
non-zero variable element has a value of uq is 0, and the probability that detS'„ = and at 
least one term in the expansion of det Sn via S„ is non-zero is 0. End of Proof of Claim 2. From 
Claims 1 and 2, it follows that if Sn is of type A„ or B„, then the probability that detM„_|_i 
undergoes the least change upon perturbation of clement rriij is equal to the probability that 
every term in the expansion of det Sn is 0. Now consider a matrix Sn of type C„. We will refer 
to the term sn • • • Snn in the expansion of det Sn via S„ as the diagonal term, and to each of 
the other terms as a non-diagonal term. Note that the diagonal term cannot be 0, since all 
the diagonal elements of Sn are fixed at 1. If xq is a number other than 1, then by reasoning 
analogous to that which was used in the proof of Claim 1, the probability that det Sn = xoisO. 
Also, by reasoning similar to that which was used in the proof of Claim 2, the probability that 
det Sn = ^ and at least one non-diagonal term in the expansion of det Sn via S„ is non-zero 
is 0. Thus if Sn is of type C„, the probability that detM„_|_i undergoes the least change upon 
perturbation of element niij is equal to the probability that every non-diagonal term in the 
expansion of det Sn is O.B We will use the term binary matrix to mean a matrix of O's and I's. 
For a matrix Sn of type A„, B„, or Cyi, let Sn — {^^i} t*e the n x n binary matrix with Sij — 1 
if Sij 7^ 0, and Sij = if Sij = 0. We will say that Sij is a variable (resp. fixed) element of Sn 
if Sij is a variable (resp. fixed) element of and that Sn is of type A„ (resp. B„, C„) if 5„ 
is of type A„ (resp. B„, C„). Every variable element of Sn has probability r of being 1, and 
probability 1 — r of being 0. Denote the permanent of any matrix A by per A. The expansion 
of the permanent of Sn via S„ is 

per5„= ^ Sa(l)lSa(2)2---S(r(n)n (6) 

ctGS„ 

Since Sn is a binary matrix, every term in this expansion is either or 1, so the value of 
perS'n is a non- negative integer. The next result follows from Proposition 1.1 together with 
techniques analogous to those used in the proof of it. 

Proposition 1.2. Let u be a real number. If Mn+i, Sn, and X are as hypothesized in 
Proposition 1.1, then the probability that det Sn = u is equal to the probability that perS'„ = u. 
Thus finding the probability that Mn+i undergoes the least change upon perturbation of element 
TUij (of Mn+i) is equivalent to the problem of finding the probability that per Sn = if Sn is 
of type An or B„, and to the problem of finding the probability that per Sn = I if Sn is of type 

Throughout section 2, we assume X to be a continuous set such that £ X. More will be said 
about this assumption in section 3. 



2. Continuous X 

Let Sn be a matrix of type A„, B„, or C„, so that some of its elements are variable, its 
other elements are fixed at 1, and each of its variable elements has probability r of being 1 
and probability 1 — r of being 0. In this section, we formulate the probability that per 5„ = u, 
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where u = if Sn is of type A„ or B„, and u = 1 if is of type C„. In what follows, we will 
use the expression per5„ = u to mean that perS'n = if Sn is of type A„ or B„, and that 
per 5„ = 1 if 5„ is of type C„. Also, we will say that a matrix 5„ of a given type (A„, B„, or 
C„) is a pertinent matrix of that type if per Sn = u- For a matrix Sn of a given type, let m 
be the number of variable elements in Sn, and let imax be the maximum number of variable 
elements with a value of 1 that a pertinent matrix of that type can have. Then the probability 
that per 5„ = u is 

^max 

^perS„=„W = E ^nii) -r'-il- rr-\ (7) 

where En{i) is the number of pertinent matrices of that type which have exactly i variable 
elements with a value of 1 — except for type B„, where En{i) is the number of pertinent 
matrices which have exactly i variable elements with a value of 1 and whose special diagonal 
element bn is in the same location as the special diagonal element of 5„. The values of m and 
imax depend on the type of matrix. Clearly, the value of m for type A„ (resp. B„, C„) is 
(resp. — n + 1, — n). To determine the value of imax for a given type, it may be easier 
to first find the minimum number jmin of variable elements with a value of that a pertinent 
matrix of that type can have, and then use the relation imax = m — jmin to compute imax- 
There exist pertinent matrices Sn of types A„ and B„ in which all n variable elements in a 
single row or column have a value of and all the other variable elements have a value of 1. 
It can easily be shown by induction on n that per 5„ > if Sn is of type A„ or B„ and has 
no row or column comprised of n variable elements with a value of 0. Thus jmin = n (and 
imax = m — n) for types A„ and B„. In a matrix of type A„, any of the n rows or n columns 
could be the single row or column comprised entirely of variable elements with a value of 0, so 
the total number of candidates for that row or column is n + n = 2n (except for n = 1, where 
that row and column coincide). Thus En{imax) = 2n if n > 1, and i?i(imax) = 1. In a matrix 
of type B„ with special diagonal element fen, there are just two candidates for the single row 
or column comprised entirely of variable elements with a value of (namely, row 1 and column 
1), the only exception being n = 1, where row 1 and column 1 coincide. Thus i?n(*max) = 2 
if n > 1, and -Ei(imax) = 1- If a matrix of type C„ is in upper-triangular form (meaning that 
every element below the main diagonal is 0) and all the elements above the main diagonal 
have a value of 1, then per5'„ = 1. Similarly, if a matrix of type C„ is in lower-triangular 
form (meaning that every element above the main diagonal is 0) and all the elements below 
the main diagonal have a value of 1, then perS'„ = 1. Such a matrix has (n^ — n)/2 variable 
elements with a value of 0, so jmin < {n'^ ^ '^)/2- In section 2.3 we show that jmin = ("^ — '^)/2 
for type C„. However, we will see that En{imax) > 2 if n > 2, so there are pertinent matrices of 
type C„ which are in neither upper- nor lower-triangular form but have imax variable elements 
with a value of 1. What remains to be done is to evaluate £'„(i) for each type of matrix, which 
entails counting matrices with per 5„ = u. 

2.1. Case (i): Type A„ 

Denote Sn by An, and En{i) by F„(i). Since m = and jmin = n, we see that imax = —n, 
so the probability that per A„ = is 

i=0 



5 



6 



G. ITO 



Thus to evaluate -Ppd- a„=o('')' have to find the integer sequence {Fn{i)}l^Q The values 
of Fn{i) for n = 1,...,4 (and i = 0, . . . , — n) are given in Table I. For n = 5, the 



sequence (from i = through i = nr 



52 - 5 = 20) is 1, 25, 300, 2300, 12650, 53010, 



174700, 458500, 956775, 1571525, 2010920, 1994200, 1534800, 923700, 439600, 166720, 50025, 
11500, 1900, 200, 10. See Sloane's sequence A088672 [5] for the first few terms of the sequence 



Table I. Values of F„(i) for n = 1, . . . , 4 (and i = 0, . 



n\i 
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9 


10 


11 


12 
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1 
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4 


4 






















3 


1 


9 


36 


78 


90 


45 


6 














4 


1 


16 


120 


560 


1796 


4080 


6496 


6976 


4860 


2128 


576 


96 


8 



{/„} = ^"(*)| (i-*'-' ^^'^ f®'*^ °f sequence 1, 9, 265, 27713, 10363661,. . .). 

Its asymptotic behavior has been proposed by Everett and Stein [2] to be 2("^~"+i). 



2.2. Case (ii): Type B„ 

Denote Sn by and En{i) by Gn{i)- Since m = — n + 1 and jmin = n, we see that 
«max = [n^ — n+1) — n = (n — 1)^, so the probability that per S„ = is 

P^erB^^r)^ J2 G^{zy{l-rf (9) 

i=0 

Thus to evaluate ^'perB„=o('')' ^® have to find the integer sequence {G„(z)}*~q ' . The 
values of G„(i) for n = 1, . . . , 4 (and i = 0, . . . ,{n — 1)^) are given in Table II. For n = 5, the 

sequence (from i = through i = (n - 1)^ = (5 - 1)^ = 16) is 1, 20, 186, 1056, 4035, 10836, 
21032, 30212, 32829, 27520, 18062, 9324, 3741, 1128, 240, 32, 2. At this time, the sequence 

{9n} = I Gn{i) > is not given in Sloane's list [5]. 



Table II. Values of Gnji) forn = 1, . . . , 4 (and ^ = 0, . . . , (n - 1)^) 



n\i 
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13 


10 
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63 


184 


315 


324 


203 


78 


18 
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2.3. Case (in): Type C„ 

Denote Sn by Cn, and by Hn{i). The probability that per Cj^ — 1 is 

^perC„ = l W = E " O"'"""' (10) 

Thus to evaluate -PperC'„=i('')' have to determine the value of imax and find the integer 
sequence {iJ„(i)}-Zo"". 

Proposition 2.1. T/ie number of matrices of type C„ wii/i permanent 1 and i variable 

elements with a value of 1 is equal to the number of acyclic digraphs with vertex set 
{f , 2, 3, . . . , n} and i edges. (The term digraph is short for directed graph; an acyclic digraph 
is a digraph that contains no cycles.) 

Proof For every matrix 5„ of type C„, let G{Sn) be the digraph with vertex set {1, 2, 3, . . . , n} 

and adjacency matrix Sn^In, where /„ is the nxn identity matrix. Then the edge set of G{Sn) 
is {{k,l) : k I, Ski = 1}, and the number of edges of G{Sn) is equal to the number of variable 
elements of Sn with a value of 1. Claim pev Sn = 1 if and only if G{Sn) is an acyclic digraph. 
Proof of Claim All the diagonal elements of Sn ~ In are 0, so G{Sn) has no loops (a loop is an 
edge that connects a vertex to itself). If G{Sn) has a cycle, then there exist j (with 2 < j < n) 
and distinct vertices fci, . . . , kj of G{Sn) such that (^1,^2), (^2, fcs), . . . , {kj-i,kj), {kj, fci) are 
edges of G(S'„). Thus Sfc^^jj 'Sfe2fe35 • • • ) ^fe^-iA-, ^ -Si^fci arc variable elements of S'„ that have a 
value of 1. If j — n, then Sk^k2 ' Sfe2fe3 ■ ' ' Sk„_ik„ ■ Sk„ki is a non-diagonal term in the expansion 
of perSn via &n which is equal to 1. If j < n, let li,l2, ■ ■ ■ , In-j be the elements of {1, 2, 3, . . . , n} 
which are not in the set {fci, . . . , kj}. Then Sfejfcs • Sfc^fea ' ' ' ■Sfe^-ifc, • Sk^ki ■ Shh ' Shh ■ ' • 
is a non-diagonal term in the expansion of pcrS'„ via 6„ which is equal to 1. In either case. 
perS'n > 1. If G{Sn) is acyclic, then the only non-zero term in the expansion of perS'„ via S„ 
is the diagonal term, so perSn = 1- End of Proof of Claim Clearly, G is a one-to-one function. 
Thus by the Claim and the definition of G, the number of matrices of type C„ with permanent 
1 and i variable elements with a value of 1 is equal to the number of acyclic digraphs with 
vertex set {1, 2, 3, . . . , n} and i edges. ■ Now suppose that a matrix Sn of type C„ has fewer 
than (n^ — n)/2 variable elements with a value of 0. Since Sn has a total of — n variable 
elements, it must have more than (n^ — '^)/2 elements with a value of 1. Thus there are k,l 
such that I < k < I < n, and (variable) elements Ski and sik of Sn that have a value of 1. But 
then the graph G{Sn) has a cycle of length 2, namely, (fc, /), (/. A:). By the Claim, perS'„ > 1. 
From this it follows that a pertinent matrix of type C„ cannot have fewer than (n^ — n) /2 
variable elements with a value of 0. As shown earlier, there exist pertinent matrices of type C„ 
with exactly (n^ —n)/2 variable elements that have a value of 0. Thus jmin = {n"^ — n)/2 (and 
imax = {n'^ — n) — [{n^ — n)/2] — (n^ — n)/2). By Proposition 2.1, the problem of evaluating 
Hn{i). the number of matrices of type C„ with permanent 1 and i variable elements with 
a value of 1, can be transformed into the problem of counting acyclic digraphs with vertex 
set {1, 2, 3, . . . , n} and i edges. We can use a result from graph theory to find Hn(i). I. M. 
Gessel [3] counted labeled acyclic digraphs according to the number of sources, sinks, and 
edges. (A source is a vertex with in-degree 0, and a sink is a vertex with out-degree 0.) We 
use that result in the proof of Proposition 2.2 below. 
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Proposition 2.2. 



Hnie) 

Yn{t) 

Z{z,t) 



1 d'Ynjt) 
e! ' df" 



OO 

E 



t=0 
Qn 



5^" \Z{-Z,t) 



2=0 



(11) 
(12) 

(13) 



n=0 



Proof For an acyclic digraph £>„ with vertex set {1, 2, . . . , n}, let e(-D„) and ,s(Z)„) be the 
numbers of edges and sources, respectively, of D„. Then define a function Ynit; a) as 

where the sum is over all acyclic digraphs £)„. (Throughout this proof, an expression of the 
form a!' is to be interpreted as 1 if both a and b arc 0. For example, is to be interpreted as 
1 if t = 0, and a° is to be interpreted as 1 if a = 0.) Let Pn(e, s) be the number of -D„ with e 
edges and s sources. Then 



Pn{e,s) = — 



t=a=0 



and we can express Yn{t; a) as 

Yn{t; a) = X! X!^"(^' ^) ■ ■ 

e 5 

Denoting Yn{t-^ 1) by Yn{t)^ we have 

Ynit) = ^ {j2Pn{e, S) j • = J2 Hn{e) ■ 

where Hn{e) is the number of -D„ with e edges. Thus 

1 rf-y„(t) 



What remains is to find a formula for Yn(t) that will enable us to evaluate Hn{e). The following 
relation is proved in [3]: 



^(1 + t)^("-^)„Qa^y„_, (t; 1) = Yn{t; a + 1) 



It is well known that 

n 



(14) 



(15) 



j=0 



.n=0 



n! ■ (1 + i)"C'2 



.n=0 



n! • (1 + t)"C'2 



^''"nl-fl 



.n=0 



n!-(l + t)"C'2 
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Setting a„, 6„, and c„ to a", Yn{t; 1) (= Yn{t)), and Yn{t;a + 1), respectively, we see that 
(14) is equivalent to 



E 



n! • (1 + t)"^^ n! . (1 + <) 



.C2 



E 



n! • (1 + f)"C'2 



(16) 



There is only one acyclic digraph with vertices (namely, the empty graph, which has no edges 
and no sources), and for n > 1 every acyclic digraph with n vertices has at least one source. 
Thus 

Yo{t; 0)=po(0,0)-i°-0° = l 
n>l r„(<;0) = ^ |^p„(e,s)-OM •t^ = 

e \s>l / 



Setting a to —1 in (16) and defining Z{z,t) as 



1 



n=0 



we find that 



71=0 



n\-{l + t)nC2 Z{-z,ty 



Thus we obtain (12). ■ Example We show the calculation of Hi{i) as an example. By (13), 



Z{-z,t) = 1 - z + 



+ ■ 



2(1 + i) 6(1 + 4)3 24(1 + i)6 120(1 + 4)10 

In the calculation of Yi{t) (see (12)), we take the fourth partial derivative of l/Z{—z,t) with 
respect to z and evaluate it at z = 0: 



54 



9zM 1 . 



1 



2(l+t) 6(T+tF 24(1+4)8 



120(1+*)"' 



The terms of Z{—z, t) that contain a factor of z*^ for some A: > 5 will vanish when the fourth 
partial derivative of l/Z{—z,t) is evaluated at z = 0, so it suffices to use the following finite 
sum, Z4{—z,t), in place of Z{—z,t) when we perform the differentiation: 



Z4(-Z,t) = 1 - Z + 



+ 



2(1 +i) 6(1 + 4)3 24(1 + 4) 



Thus 



F4(4) = (1+4)^ 



= + 



r M 



dz'^ \Z{-z,t) 
( 



z=0 



dz^ 



1-Z+- 



+ 



2(1 + 4) 6(1 + 4)3^24(1 + 4)6/ 
= 1 + 124 + 604^ + 15243 + 1864^ + 1084^ + 244^ 



z=0 
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That is, (i?4(0), H4{1), Hi{2), Hi{'i), H^A), H^B), ^4(6)) = (1, 12, 60, 152, 186, 108, 24). 
Thus 

^perC.=lW=E^^4(iK(l-r)12- 

1=0 

= 1 • (1 - r)i2 + 12 • ri(l - r)" + 60 • r'^^l - r)^° + 152 ■ r^{l - rf 
+186 • r^(l - rf + 108 • r^(l - rf + 24 • r^(l - rf 

■ The values of Hn{i) ioi n — 1, . . . , 5 (and i = 0, . . . , (n^ — n)/2) are given in Table III. See 
Sloane's sequence A003024 [5] for the first few terms of the sequence = {^^i?„(i)} (i.e., 

the first few terms of the sequence 1, 3, 25, 543, 29281,. . .). As shown in the table, there are 



Tabic III. Values of H„{i),n = 1, . . . , 5 (and ^ = 0, . . . , (n^ - n)/2) 



n\i 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


1 






















2 


1 


2 




















3 


1 


6 


12 


6 
















4 


1 


12 


60 


152 


186 


108 


24 










5 


1 


20 


180 


940 


3050 


6180 


7960 


6540 


3330 


960 


120 



six matrices of type C3 with permanent 1 and imax (= 3) variable elements that have a value 
of 1. Of those six matrices, there are four that are in neither upper- nor lower-triangular form: 

101\ /100\ /111\ /110\ 
1 1 1.1 1 1, 1 , 1 

001/ \ioi/ \oiiy \iiiy 

2.4- Remarks 

Graphs of the functions ^'pcj. a5=o(^)' -^pori3r,=o('')' ^pcrC5=i('') presented in 

Fig. 1. If we were to start with a matrix of type A„ (case (i)) and replace all but one of 
the diagonal elements with a fixed element 1, we would obtain a matrix of type B„ (case 
(ii)). If wc in turn replaced the special diagonal clement in the latter matrix (the variable 
element 611) with a fixed element 1, we would obtain a matrix of type C„ (case (iii)). 
Therefore, we might expect to find that there exist analytic formulas for ^'per b„=o('') terms 
-Ppcr A„=o('')' ^^"^ ^pcrC'„=o(^) tcrms of Ppa b^=q{''^) ■ At thc vcry least, wc might expect 
that thc relation Pp^r A„=oi^) > -PperB„=o(^) > -PperC'„=i(^) would hold and that thc value 
of f'pcr A„=o(^) - Ppcrsl^oir) would bekrgcr than the value of fpcrB„=o(^) -^pcrC'„=i('')- As 
can be seen in thc figure, such a simple relation holds for n = 5, but only for values of r larger 
than approximately 0.18. There is another relationship between matrices of type B„ and those 
of type C„+i. It can be shown that if row i and column j are deleted from a matrix of type 
C„+i, thc resulting n x n matrix can be converted to a matrix of type B„ by interchange 
of at most one pair of rows or one pair of columns. It is well known that the absolute value 
of the determinant of a matrix is invariant under interchange of any two rows or any two 
columns. Thus the cofactor of any element of a matrix of type C„+i is equal in absolute value 
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P [r. 5] 



Figure 1. (i) PperA,=o(0, (") PperB,=oir), (i") P^erC,=iir) 

to the determinant of a matrix of type B„. Though we obtained a closed formula for the nth 
term of the sequence {hn} = {'^Hn{i)} (where hn is the total number of matrices of type 

i 

C„ with permanent 1), we did not obtain closed formulas for the nth term of the sequences 
{/„} = {'^Fn{i)} and {gn} = (where /„ is the total number of matrices of type 

i i 

A„ with permanent 0, and g„ is the number of matrices of type B„ with permanent and a 
fixed location of the special diagonal element 6ii). We found that it took considerably more 
time to compute {/„} and {gn} than to compute {/i„}. It is possible that the problems of 
computing {/„}, {gn}, and {hn} arc related to some ^P- complete problem. ^P-completc is 
a class of counting problems in complexity theory. The problem of computing the permanent 
of a binary matrix is ij^P- complete, as reported by Valiant [6] and by Ben-Dor and Halevi [1]. 
The determinant of a given matrix can be computed in cubic time by Gaussian elimination, 
but that approach cannot be used to compute the permanent. The key difference between 
determinants and permanents is that the relation detAB = detA ■ detB always holds but the 
relation peiAB = per A ■ peiB does not hold in general [4]. 



3. Continuous X vs. discrete X 

Let m(A„, X) be the real number u of least absolute value such that a matrix of type A„ can 
have determinant u with non-zero probability if its variable elements are chosen from X. If 
such u exists and X is discrete, let f2(A„, X) be the set of all matrices Sn of type A„ such that 
det Sn = u{An,X) with non-zero probability. If such u exists and X is continuous, define the 
following equivalence relation on the set of matrices A„ of type A„ such that det An = w(A„, X) 
with non-zero probability: Matrices M and N are equivalent if the locations of the variable 
elements that have a non-zero value are identical in M and N. Then let f2(A„,X) be a set 
that contains one representative of each equivalence class. For every matrix Sn of type A„, let 
Sn be the corresponding binary matrix of type A„, i.e., the binary matrix with element S,j 
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defined by 

r 0, ifsi, = o 

\ 1, if Si,- 7^0 

Let every variable element of a binary matrix Sn of type A„ have probability r of being 1, 

and probability 1 — r of being 0. Then let w(A„,X) be the real number u of least absolute 
value such that a matrix of type A„ can have determinant ii with non-zero probability if the 
variable elements of matrices of type A„ are chosen from X, and let n(A„, 3£) be the set of all 
matrices Sn of type A„ such that det S*,! = u(A„,X) with non-zero probability. In addition, 
for every i < imax, let r2j(A„,X) be the subset of f2(A„,X) that consists of matrices of type 
An that have i variable elements with a non-zero value, and let f2i(A„,3£) be the set of all 
matrices in r2(A„,X) that have i variable elements with a value of 1. (imax is the maximum 
number of variable elements with a value of 1 that a matrix of type A„ which has determinant 
u{An,3i) with non-zero probability can have.) For type B„, define the following in a similar 
manner: m(B„,X), fi(B„,X), 'u(B„,X), f2(B„,X), and, for every i < imax, fii(B„,X) and 
rii(B„,X). Then do likewise for type C„. For example, u(C2, [0,2]) = 1. and 

o,c.io..i,={(i ;)}. m 

where x and y are elements of 3£ — {0}. The probability that a matrix C2 of type C2 is 
(^0 1 ^ ~ ^^^^^ probability that C12 = is 1 — r, as is the probability that 

/I a; \ 

C21 = 0. The probability that C2 is in the equivalence class of I 1 is r • (1 — r), because 



^0 1 ^ 

the probability that C12 7^ is r and the probability that C21 = is 1 — r. Similarly, the 
probability that C2 is in the equivalence class oi 1 ^ ~ '')• example of the 

discrete case, let 3C = {0, 1/2, 2}, and let p(l/2) = r/2 = p(2). Then a matrix C2 of type C2 

can have determinant with non-zero probability (i.e., m(C2, {0, 1/2, 2}) — 0), because the 
probability that the term C12C21 in the expansion of detC2 via 62 is equal to 1 is non-zero. 
Note that ^ 

J^(C2,{0,1/2,2})=| ( 1 , ), ( ^ 2)1. (19) 



The probability that C2 = ( \ ? ) is {r/2)^, as is the probability that C2 = ( ^ 2 ) • 
Thus 



and 



O(C2,{0,l/2,2})= j ( J J 1 ^ (21) 



Partitioning the sets f7(C2, [0, 2]), 17(C2, {0, 1/2, 2}), ii(C2, [0, 2]), and $1(C2, {0, 1/2, 2}) 
according to the number i of variable elements which arc non-zero, we find that 

f^o(C2,[0,2]) = {( J Ml, f]i(C2,[0,2]) = {( J J)} (22) 
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Oo(C2,[0,2]) = |(^ J J )|, Oi(C2,[0,2]) = |(^ J \ J)} (23) 

and 

Oi(C2,{O,l/2,2}) = 0, z = 0,l, f22(C2,{0,l/2,2})= If } 2^1 (24) 



fji(C2, {0, 1/2, 2}) = 0, i = 0, 1, ^2(62, {0, 1/2, 2}) = I ( J J j I (25) 

For sets S, T, we will use 5 ~ T to denote that there is a one-to-one correspondence between 
S and T. 

Proposition 3.1. If X is continuous, then 
1. ■u(A„,X) = w(A„_,X) 

3. the following are equivalent for a matrix Sn of type A„; 

(a) detS'n = m(A„,X) with non-zero probability 

(h) perSn — u{AmX) with non-zero probability 

(c) det^n = u{An,X) with non-zero probability 

(d) per^ra = {((A„,X) with non-zero probability 

The counterparts of all of these statements hold for types B„ and C„ as well. IfXis discrete, 

then statement 1 holds for types A„ and B„, but it does not hold in general for type C„; also, 
neither statement 2 nor statement 3 holds in general (for any of the three types). 

Proof We have already established this for the case where X is continuous: For types A„ 
and B„, we showed that the least value of u such that det5„ = u with non-zero probability 
is 0, and that dct S'„ = with non-zero probability if and only if per = with non-zero 
probability. For type C„, we showed that the least value of u such that det Sn = u with non- 
zero probability is 1, and that det 5n = 1 with non-zero probability if and only if per Sn = 1. 
The truth of statement 2 is obvious. For type A„ or B„, the equivalence of 3(a) and 3(b) 
follows from our proof that det Sn = with non-zero probability if and only if every term in 
the expansion of det 5„ via (3„ is 0, which is true if and only if every term in the expansion of 
pel Sn via 6„ is 0. For type C„, the equivalence of 3(a) and 3(b) follows from our proof that 
det S'n = 1 with non-zero probability if and only if every non-diagonal term in the expansion 
of det Sn via (S„ is 0, which is true if and only if every non-diagonal term in the expansion 
of per S'„ via (5„ is 0. If X is discrete, then for types A„ and B„, the least value of u such 
that det Sn = u with non-zero probability is 0, as is the least value of u for which det Sn = u 
with non-zero probability. Both of these facts are witnessed by matrices that have a single 
row or column comprised entirely of variable elements with a value; of 0. For a counterexample 
to statement 1 in type C„, let X = {0, 1/2}. There are only four matrices of type C2 with 
variable elements chosen from X: 

(if) (n) (n 
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Thus m(C2,3£) = 3/4, as witnessed by ^2 = ^ i f ^' ^^i^h has determinant 3/4 with 

non-zero probability. Tlic associated binary matrix is ^2 = ^ | | ^ ' ^^i^^li ^'^^ determinant 

with non-zero probabihty, so {((C2,X) = 0. For a counterexample to statement 2 in the 

1 




discrete case, let 3C = {0, 1/2, 1, 2}. If = | q ^ q | , then Ss can be viewed as a matrix 

2 11, 

/ 1 1 

of type A3, type B3, or type C3. The associated binary matrix is ^3 = 1 | . Since 

. V 1 1 1 

{0, 1} C X, the 3x3 matrix T3 whose elements are equal to those of S3 can also be viewed as 
a matrix of type A3, type B3, or type C3. Regardless of type, matrices S3, T3, and 6*3 all have 
determinant with non-zero probability, so u and u are (even for type C3). Thus ^3 and 
T3 serve as witnesses to the failure of statement 2: The fact that S3 is the associated binary 
matrix for both ^3 and T3 means that fie (A3, X) 9^ J76(A3,X), 1^4(63, X) 9^ Qa{B3,3l), 
and $73(03, X) 9^ $73(03, X). Also, per 53 = 2 = per 53. Thus for all three types of matrices, 
5*3 also serves as a coimterexample to statement 3, though not as a counterexample to the 
equivalence of 3(a) and 3(c), and not as a counterexample to the equivalence of 3(b) and 3(d). 
For a counterexample to the equivalence of 3(a) and 3(c) (for all three types of matrices) in 
the discrete case, let X = {0, 1, 2}. Then u = = u for types A3, B3, and C3, as witnessed by 
/I 1 \ 

the 3x3 matrix M = 1 1 I , which has determinant with non-zero probability and 

\0 1 1] 

can be viewed as a matrix of type A3, type B3, or type C3. Moreover, since {0, 1} C X, the 
elements of the associated binary matrix of M are identical to those of M itself, so det M = 

/I 1\ 

with non-zero probability. If ^3 = I 1 1 , then S3 can also be viewed as a matrix of 

V 2 1 1 / 

type A3, type B3, or type C3. The associated binary matrix is >S'3 = I 1 1 . Thus 

V 1 1 1 / 

dot 5*3 = and dot S3 — 1, which shows that 3(a) and 3(c) arc not equivalent in general. Hence 
3(b) and 3(d) are the only two items in statement 3 whose equivalence we have neither proved 
nor refuted (in the case where X is discrete). ■ 

Proposition 3.2. Let X^™*) and X^^'") be a continuous set and a discrete set, respectively, 
such that e X^^^^) C X^™*). Then f2(A„, X^^^^^)) □ Q(A„,X('^"*)) and n(B„,X(d^^)) □ 
r2(B„, X^'^"*-'), where C and C denote set inclusion for sets of numbers and sets of matrices, 
respectively. However, there exist X^'^"*' and X^*^"*) such that e X^**"*) C X^™'^^ but 
0(C„,X(d-))2 0(C„,X(-*)). 

Proof The first statement follows from the fact that every element of f2(A„, X^'^'^*)) IS a 
matrix with at least one row or column comprised entirely of variable elements with a value of 
0, as is every element of ii(B„, X(™*)). For example, let X^'^"*) = [0, 2] and X(<i'^) = {0, 1/2, 2}. 
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Then u(A2,X(™*)) = = u(A2,X(<^'")) (therefore m(A2,X(™*)) = = u(A2, so 



n(A2,x('^"*)) = I 



00\/10\/0l\/00\/00 

oj'l^o ) ' y 1 oj'lo 1 

1 l \ /O 0\ / I OWO 1 

oi'li ij'li oilo 1 



f](A2,X('"^)) = 51(A2,X('="^))|J I 



1 1 
1 1 



(26) 



(27) 



The example given in (20) and (21) illustrates the truth of the second statement. ■ Wc 
can exploit the difference between the continuous and discrete cases to create an interesting 
example in which the real number u of least absolute value such that a matrix of typo C2 can 
have determinant u with non-zero probability in the continuous case differs from its counterpart 
in the discrete case. Let X^™*) = [0, 1] and X(<i'^) = {0, 1}. Then u(C2,X(™*)) = 1 (therefore 
{t(C2,3e(^"*)) = 1), and 



n(C2,3e('^"*') : 
n(C2,3e(™*)) 



1 \ / I X \ I 

I I )'\y 1 

1 o\ fi i\ fi 
1 / ' I 1 / ' I 1 1 



where x and y are elements of X^'^"*) -{0}. Also, m(C2, X^^^'")) = (therefore u(C2, X^^^'")) = 0), 
and 



f^(C2,X(d-)) = [^ \ \ )| =i^(C2,X(<^-)) (28) 



Thus 

fi(C2, x(™*)) n o(C2, x(<^'")) = 0. 

The probability that detC2 = 1 (perC2 = 1) in the continuous case is 

^iTc.=iM - - rf + 2r(l ^ r) = 1 - r^ 

while in the discrete case the probability that det C2 = is 

Therefore, 

P^^k=^ir)+P^:l^,{r) = l (29) 
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